Magnetic winding and method of making same

ABSTRACT

The present invention provides an improved magnetic winding and method of calculating desired winding parameters (winding layer thickness, number of winding layers and number of turns per winding layer) for a winding in a magnetic component. The invention may be applied to general boundary conditions in a magnetic winding or component and considers relative phase displacement for sinusoidal and nonsinusoidal winding currents. Ratios of magnetic surface field intensities at corresponding inner and outer boundaries of one or more winding layer(s) are calculated, and considered with relative phase displacement to select magnetic winding configurations having desired or optimal power dissipation. In certain aspects, a normalized loss function f(H,R,B,Φ) is utilized to determine a preferred construction among a plurality of iteratively generated selections.

TECHNICAL FIELD

The present invention relates generally to magnetic components and, more specifically, to an improved method of calculating or selecting desired magnetic winding parameters including winding layer thickness, number of winding layers and number of turns per winding layer, and an improved wound magnetic component such as a coil, inductor, transformer or motor having reduced power dissipation.

BACKGROUND OF THE INVENTION

Magnetic and inductive components and devices such as coils, inductors, transformers, motors and others (referred to herein as magnetic components) include a winding of one or more of a variety of conductors. Such magnetic components utilize a variety of conductor types including round, square, or rectangular wire; thin conductor strips or foil; multiple wires twisted together or wound in parallel; various Litz or woven magnet wire braids (to increase uniformity of current sharing between elementary conductors); and various combinations of such conductor types.

Conventional solenoidal magnetic components comprise winding layers which form cylinders around a core. In contrast, planar magnetic components utilize conductors and combinations of conductors in an annular configuration. Round, square or rectangular wire, for example, may be wound in annular configurations. In addition, thin conductor strips or foil may be implemented in annular configurations using printed circuit boards, flex circuits, or discrete conductors fabricated from sheet stock, for example, and in other configurations known to those skilled in the art. In comparison to solenoidal magnetic components, the thickness of winding layers and the number of turns per layer in planar magnetic components, such as those on printed circuit boards, may be varied easily and inexpensively.

FIG. 1 illustrates a cross-section of a winding region or portion in a typical magnetic component having n layers. In FIG. 1, the winding length of each layer is designated l. The thickness of each winding layer is designated T₁ through T_(n) and each layer has N_(i) turns, designated N₁ through N_(n). The magnetic surface field intensities at the inner and outer boundaries of the ith winding layer are designated as H_(i-1) and H_(i) respectively. The current in the ith winding is designated I_(i) and points out of the plane of the paper. When the winding length l is much greater than the winding layer thickness, the magnetic field distribution is largely parallel to the plane of the conductor in each winding layer. The magnitude ratio of peak magnetic surface field intensities for each conductor layer is defined as follows: R_(n)=H_(n)/H_(n-1) for each of n layers. Phase shift or phase displacement of magnetic surface field intensities for each conductor layer is defined as: Φ_(n)=φ_(n)−φ_(n-1). The turns in FIG. 1 are illustrated by way of vertical lines in each of the winding layers.

Unless otherwise stated, the following additional definitions with implied units are used herein:

-   -   H: Magnetic Field Intensity: units of

$\frac{\text{Ampere-Turn}}{\text{Meter}}$

-   -   ρ: Resistivity: units of ohm—meter     -   μ₀:

$\text{Permeability~~constant} = {4\pi \times 10^{- 7}\frac{Henry}{\text{Meter}}}$

-   -   f: Excitation frequency: units of Hertz     -   δ:

$\text{Skin~~Depth} = \sqrt{\frac{\rho}{{\pi\mu}_{0}f}}$

Magnetic components always incur some power dissipation in the winding(s) and core, which decreases efficiency and increases temperature. It is generally known that alternating current (AC) conduction generates eddy currents within the conductors of magnetic components. Such eddy currents are significant at high frequencies and/or for large conductor thicknesses. These eddy currents do not contribute to the macroscopic current of the device, but produce a field which tends to cancel the external magnetic field produced by the AC current. However, the resultant power dissipation, or loss, and energy storage associated with such eddy currents can have a significant impact on the performance of a magnetic component in an electrical circuit. In particular, dissipation from eddy currents can markedly reduce the electrical efficiency of a system and increase the temperature rise of the component. This is due to the well known skin and proximity effects.

Skin effect is the tendency of the current density in a wire to increase at and near the surface of the wire. In other words, skin effect is the tendency of current in a conductor to flow more toward the surface of the conductor as frequency is increased. Current density decays exponentially inside the conductor, reaching a value at the skin depth (δ) of 1/e times the current density at the surface.

Proximity effect occurs when one conductor is placed in an external field generated by one or more other conductors in close proximity. In that case, eddy currents are induced in the conductor which oppose the penetration of the external field. The two eddy current effects occur simultaneously in a conductor carrying an AC current when the conductor is exposed to an external magnetic field. Such eddy currents cause power dissipation in the windings of magnetic components which increase with frequency and/or at large conductor thicknesses.

Heretofore, designers of wound magnetic components have utilized analytical methods to limit or reduce power dissipation based upon mathematical derivations by P. L. Dowell in 1966 (P. L. Dowell, “Effects of Eddy Currents in Transformer Windings,” Proceedings of the IEEE, Vol. 113, No. 8, August 1966 [incorporated herein by reference]). P. L. Dowell, as with most prior art, assumes a constant layer thickness or height and does not consider current phase displacement.

High frequency analysis of coil regions has been analyzed using a classical equivalent foil representation of a winding layer. This approach facilitates an understanding of physics and determination of conductor boundary conditions. However, this method neglects stray field effects at edges and other asymmetries. Indeed, stray effects can also arise from unpredictable manufacturing variables such as insulation build up and winding terminations which can cause irregular conductor geometries.

Increased computing power has facilitated iterative approaches to determine winding configurations that yield acceptable dissipation. For example, Finite Element Analysis has improved mathematical consideration of asymmetries and specific device geometries. Finite Element Analysis is frequently used to determine whether a specific component design is acceptable. It is less valuable, however, in generating or suggesting all potential design parameters and determining an optimal or desired solution. Finite Element Analysis software may examine the impact of various configuration parameters such as core type, conductor type and size and the configuration of terms without the need to build and test a physical device.

Prior to the present invention, the minimum loss configurations for the individual winding layers of a magnetic component having more than one layer was not analytically derived for the general case. M. P. Perry, “Multiple Layer Series Connected Winding Design for Minimal Losses,” IEEE Transactions on Power Apparatus and Systems, Vol. PAS-98, No. 1, January/February 1979, pp. 116-123, discloses an analysis for minimized power dissipation by choosing specific radial thicknesses for each winding layer in a magnetic component. M. P. Perry's analysis is based on general field solutions for current density distributions in winding layers of an infinitely long, cylindrical current sheet. Also, Perry's analysis assumes a fixed number of turns per layer and zero phase displacement. The analysis, therefore, has limited applicability. In addition, the stated twelve percent reduction in power dissipation in the M. P. Perry paper has subsequently been considered too small a benefit when increased manufacturing costs are considered.

Since prior art methods have focused on the equivalent AC resistance of a complete winding portion or region, the optimization or minimization of loss in discrete winding layers within a magnetic component has not been implemented. As a result, configurations of minimum winding dissipation have been elusive and magnetic components have been less efficient and larger or hotter in comparison to results for an optimized configuration.

U.S. Pat. Nos. 6,455,971 and 6,758,430 (Palma et al.) disclose a non-random winding technique to reduce proximity losses in motors and other electric machines. U.S. Pat. No. 6,617,665 (Farcy et al.) discloses an optimized width for inductive windings on an integrated circuit, the width being twice the “skin thickness” corresponding to the maximum frequency of a high frequency current running through the winding. U.S. Pat. No. 6,650,217 (Wolf et al.) discloses a low profile planar magnetic component having a stacked winding configuration and specifying a minimum distance between winding layers and an air gap. U.S. Pat. No. 6,661,326 (Yeh et al.) discloses a wire-winding structure and method to improve transformer power which consists of a method of winding wire from a pin on a bobbin around a plurality of slots. U.S. Pat. No. 6,536,701 (Fulton et al.) discloses the use of an improved former for winding electrical coils. None of these patents discloses a device or method having improved winding parameters as in the present invention.

Generally, the prior art assumes uniform conductor thickness; approximates dissipation in terms of an equivalent AC/DC resistance ratio for the entire winding portion; and suggests that the best way to reduce winding eddy current loss in a magnetic component is to reduce the number of layers. In addition, none of the prior art accounts for phase displacement in determining preferred winding parameters.

There is a need, therefore, for an improved magnetic component having desired or optimal winding parameters including winding layer thickness, number of winding layers and considering variable turns per winding layer. There is also a need for a cost-effective method of designing and/or manufacturing such magnetic components to reduce power dissipation. Since multiple transformer secondaries, for example, can have loads with unequal power factors, or differing nonlinear loads (e.g. independent secondary regulators), and since significant magnetization currents can occur in primary windings, the harmonic components of winding currents can have significant relative phase displacement. It is desirable, therefore, to provide a single method of designing or calculating winding parameters which may be applied to general boundary conditions, including consideration for relative phase displacement of winding currents.

BRIEF SUMMARY OF THE INVENTION

With reference to the corresponding steps, parts, portions or surfaces of the disclosed embodiment, merely for purposes of illustration and not by way of limitation, the present invention provides an improved method of determining, calculating and/or designing winding parameters for use in a magnetic component, such as the desired thickness of a winding layer, the desired number of winding layers and the desired number of turns in each winding layer. The novel method may be applied to a single winding layer in one-layer or multi-layer components, or to any number, or all, of the winding layers in such components. The method can also be applied to separate and distinct windings that are placed on a given layer.

The selection of a desired winding parameter, such as winding layer thickness, turns per layer, or number of layers, in accordance with this invention, may result in optimal power dissipation and/or desired power dissipation given manufacturing constraints such as size, cost, minimum/maximum leakage inductance or minimum/maximum capacitance, or other factors. By improving and/or optimizing conductor thicknesses in all or selected layers of a winding, for example, winding dissipation is significantly reduced, yielding higher efficiency and a corresponding opportunity for size or temperature rise reduction. As a result, coils, inductors, transformers, motors and various other magnetic components, and the systems which incorporate them can be made smaller and more efficient in accordance with the present invention, without undue complexity or expensive component materials.

As described above, FIG. 1 is a cross sectional view of a winding region of a magnetic component. AC current within the conductors creates gradients of AC magnetic field within the winding region. These AC magnetic field gradients apply corresponding magnetic field boundary conditions to the conductor elements which induce eddy currents which can significantly increase dissipation. In certain aspects of the invention, desired or optimal conductor geometry is determined as a function of ratio(s) of magnetic field boundary conditions such as peak magnetic surface field intensities. In other aspects, conductor/winding geometry is determined as a function of relative phase displacement between the magnetic fields at the winding layer boundaries. In one aspect, the present invention utilizes a determination of optimum or desired conductor thickness which minimizes winding layer dissipation. In another aspect, this invention specifies the optimum conductor thickness for each layer of the inductive device. Optimum conductor thickness can be implemented, however, on each layer or in selected layers based upon preferences (or limitations) in component material and manufacturing processes.

One aspect of the invention provides the following method of calculating the desired thickness of a winding layer for use in a magnetic component: determining the magnetic field intensity at the inner boundary or surface of a first winding layer; determining the magnetic field intensity at the outer boundary of the first winding layer; calculating a ratio of the magnetic field intensities; and calculating the desired thickness of the first winding layer as a function of this ratio. In certain aspects of the invention, the calculation of a desired thickness minimizes the normalized power loss function, f(H,R,B,Φ), defined in FIG. 3, wherein B is a ratio of the desired thickness to the skin depth (δ) of a winding layer, R is the ratio of the magnitudes of opposing peak magnetic surface field intensities, and Φ is the relative phase displacement between the phases of the magnetic fields at the inner and corresponding outer boundaries or surfaces of a given winding layer. FIG. 2 illustrates such opposing magnetic surface field intensities for an example winding layer (10). The normalized loss function was generated by finding the general solution of the magnetic diffusion equation. A conductor layer is represented as a semi-infinite current sheet having magnetic field boundary conditions that are parallel to the conductor surface. This assumption is valid at least when winding width, l in FIG. 1, is much larger than winding layer thickness. The application of Maxwell's equations generates the diffusion equation which has been classically solved using imaginary representations of magnetic field intensity and current density. True results are given by the real components of the imaginary solutions. The present invention is based in part upon a further derivation of P. L. Dowell's loss expressions to enable effective determination of the optimization function for desired reduced loss configuration(s) for zero phase disparity. Assumed boundary conditions were expanded to encompass phase disparity of magnetic fields and new loss expression(s) and their corresponding optimization functions derived, whose solutions define desired reduced or minimized loss configurations.

In certain aspects of the invention, when phase displacement is zero and winding current is sinusoidal (and R>1), the desired thickness of a winding layer may be calculated in accordance with the expression, cos h(B)=R*cos(B). In another aspect, the calculation of desired thickness comprises: computing power dissipation for a plurality of predetermined thicknesses (e.g. for a selected range, thickness may be varied by a selected increment), which may be limited by manufacturing constraints or otherwise; and selecting a desired thickness having a desired power dissipation. In other aspects of the invention, the method contemplates that various winding parameters be varied in the same manner. In other aspects of the invention, the calculation of a desired thickness of a winding layer comprises: plotting the power loss function, f(H,R,B,Φ), for a plurality of predetermined thicknesses; and selecting a desired thickness having a desired power dissipation. Such predetermined thicknesses may be generated by a computer or determined manually, for example.

Another aspect of the invention provides the following: determining the magnetic field intensity at the inner and outer boundaries of a second winding layer; calculating a ratio of those magnetic field intensities; and calculating the desired thickness of the second winding layer as a function of that ratio. The same determinations and calculations may be made for any number of winding layers and for one, some or all layers in a given component. The desired thickness may minimize the normalized loss function, f(H,R,B,Φ). In another aspect, the desired thickness of a second winding layer is calculated in accordance with the expression: cos h(B)=R*cos(B) as defined above. As in a single layer computation, power dissipation may be computed for a plurality of predetermined thicknesses; and a desired thickness may be selected which has a desired power dissipation. Similarly, another aspect of the invention provides a method by which the power loss function, f(H,R,B,Φ), is plotted for a plurality of predetermined thicknesses of a second (or greater) layer; and a desired thickness is selected which has a desired power dissipation.

In other aspects, the desired thickness calculated or selected pursuant to the invention is the optimal thickness of the winding layer(s). In one aspect, the method of this invention includes: manufacturing a magnetic component having such desired thickness. Another aspect provides a magnetic component having a magnetic winding or windings made according to the method(s) described above. In certain aspects, for sinusoidal winding current corresponding to a situation in which the ratio R=0, the desired thickness of the winding layer(s) may be calculated in accordance with the expression: B=π/2. When R=−1 and phase displacement is zero, the desired thickness(es) may be calculated in accordance with the expression: B=π.

In one aspect of the invention, the desired thickness of the first and/or second winding layer in a winding in a magnetic component is calculated as a function of the relative phase displacement between the phases of the magnetic fields at the inner and corresponding outer boundary of individual winding layers. One aspect of the invention provides the following method of calculating the desired thickness of a winding layer for use in a magnetic component: determining the magnetic field intensity at the inner surface or boundary of a first winding layer; determining the phase of the magnetic field at the same inner boundary; determining the magnetic field intensity at the outer boundary or surface of the first winding layer; determining the phase of the magnetic field at the same outer boundary; calculating a ratio of the magnetic field intensities determined and the relative phase displacement between the phases of the magnetic fields determined; and calculating the desired thickness of the first winding layer as a function of the ratio and the relative phase displacement. In one aspect, for sinusoidal winding current and cos(Φ)>0, the desired thickness of a winding layer may be calculated in accordance with the expression: [R ²+1]/R=[cos(B)/cos h(B)+cos h(B)/cos(B)]*cos(Φ)

For the general case with sinusoidal winding current, the desired thickness of a winding layer may be calculated in accordance with the expression: R cos(Φ) sin(B)[sin h(B)[4 cos²(B)+1]+sin h(3B)]−[R ²+1]sin(2B)sin h(2B)=0

Each of the foregoing steps may be undertaken, in accordance with the present invention, for a second winding layer or any number of winding layers in the winding of the magnetic component. In certain aspects, power dissipation is computed in accordance with the power loss function, and the power loss function may be plotted for a plurality of potential winding thicknesses, and a desired thickness selected based on the plot. The conductor which comprises the winding layer in each aspect of the invention may be round wire, square wire, foil sheet, conductive tape or multiple-strand wire. The thickness of round wires can be taken as

$\sqrt{\frac{\pi}{4}}*$ diameter so that equivalent cross-sectional areas are established for the purpose of analysis. The invention also includes the step of manufacturing a magnetic component (having a magnetic winding) in accordance with the foregoing steps. A magnetic component which is manufactured in accordance with the foregoing steps/methods is also part of the present invention.

For nonsinusoidal winding current (which may include a DC component), the invention provides for the determination, through Fourier decomposition or another method known to those skilled in the art, of one or more of the harmonic components of the winding current. For one or more of the harmonic components, calculations may be made of H, the ratio R and the relative phase displacement for the inner and outer boundaries of each considered winding layer; and a desired thickness determined from such calculations. Dissipation may be determined as a function of H, R, B and Φ in accordance with the power loss function. In another aspect of the invention, potential winding layer thicknesses are iteratively generated (by computer or otherwise), and a desired thickness may be selected which has a desired harmonic dissipation. In certain aspects, similar to the sinusoidal case, the power loss function is plotted and a desired thickness is selected by computer or by inspection. The foregoing steps may be followed for one, some or all considered winding layers.

The present invention also provides for a determination of the following additional winding parameters: number of winding layers and number of turns per layer. In one aspect, the invention provides the following method for calculating such parameters: determining the harmonic components of the winding current; iteratively generating a plurality of combinations of one or more of the winding parameters (including winding layer thickness); for one or more of the harmonic components, computing power dissipation for a plurality of the generated combinations as a function of the relative phase displacement and/or as a function of the ratio of magnetic surface field intensities between the inner boundary and the outer boundary of each of the iteratively generated winding layers; comparing the resulting power dissipations; and selecting a desired combination of winding parameters having desired power dissipation. The combinations selected may be limited by manufacturing constraints or specification constraints. As in other aspects of the invention, power dissipation may be determined in accordance with the power loss function, f(H,R,B,Φ), and the comparison of resulting dissipations done by plotting the function. These steps may be undertaken for one, some or all considered winding parameters.

In another aspect, the present invention provides for a determination of winding thickness, number of winding layers and number of turns per layer as follows: varying one or more of the winding parameters; determining boundary conditions at the inner boundary and the outer boundary of one or more winding layers for a plurality of combinations of the varied winding parameters; comparing the boundary conditions at the inner and outer boundaries; and computing power dissipation as a function of the relationships between the compared boundary conditions. Each of the foregoing steps and methods may be implemented through a computer readable medium having computer executable instructions. Such computer readable mediums include without limitation floppy and hard disks, CD-ROM, flash ROM, nonvolatile ROM, DVD, and RAM, for example.

The general object of the invention is to provide a magnetic component (and magnetic winding) having improved, desired or optimal power dissipation. This improvement can be made in consideration of imposed specification or manufacturing constraints. It is a further object of the present invention to provide a method for determining winding parameters (winding layer thickness, number of layers, number of turns per layer) for a wound magnetic component having a desired or optimal winding loss. This desired loss can be achieved in consideration of imposed specification or manufacturing constraints. Another object is to provide a method of winding design which is applicable to any type of component geometry, such as planar transformers, solenoidal transformers, and motor geometries.

Still another object is to provide a method which may be applied to general boundary conditions, including consideration for relative phase displacement of winding currents. These and other objects and advantages will become apparent from the foregoing and ongoing written specification, the accompanying drawings and the appended claims.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a cross-sectional view of a winding region of a magnetic component.

FIG. 2 is a cross-sectional view of a winding layer of a magnetic component.

FIG. 3 is an expression for normalized power loss.

FIG. 4 is a chart illustrating the application of the invention to a magnetic component with three windings.

FIG. 5 is a chart illustrating the application of another embodiment of the invention to a magnetic component with two windings.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

At the outset, it should be clearly understood that like reference numerals are intended to identify the same structural elements, portions or surfaces consistently throughout the several drawing figures, as such elements, portions or surfaces may be further described or explained by the entire written specification, of which this detailed description is an integral part. Unless otherwise indicated, the drawings are intended to be read (e.g., arrangement of parts, proportion, degree, etc.) together with the specification, and are to be considered a portion of the entire written description of this invention. As used in the following description, the terms “horizontal”, “vertical”, “inner”, “outer”, “up” and “down”, as well as adjectival and adverbial derivatives thereof (e.g., “horizontally”, “rightwardly”, “upwardly”, etc.), or similar terms, simply refer to the orientation of the illustrated structure as the particular drawing figure faces the reader. Similarly, the terms “inwardly” and “outwardly” generally refer to the orientation of a surface relative to its axis of elongation, or axis of rotation, as appropriate.

In a first preferred embodiment of the present invention, corresponding to AC sinusoidal excitation of a magnetic component having a fixed number of winding layers and a fixed number of turns per winding layer, one or more of the following steps are performed to calculate or determine a desired winding layer thickness in one or more magnetic windings:

-   -   a. Identify the current magnitude and phase for each winding in         the magnetic component.     -   b. Using Ampere's Law, calculate the boundary conditions of peak         magnetic surface field intensity and phase at each winding layer         boundary.     -   c. For each winding layer, determine the ratio of opposing peak         magnetic surface field magnitudes (R_(n)) and the relative phase         displacement, Φ_(n)=φ_(n)−φ_(n-1).     -   d. For each winding layer, determine the value of B         (B_(optimum)) which minimizes the loss function, f(H,R,B,Φ), set         forth in FIG. 3. After applying the loss function, optimum or         desired conductor thickness is then equal to B_(optimum)*δ (skin         depth).     -   e. Use one or more of the following methods to find the minimum         or desired values of the loss function:         -   i. Plot f(H,R,B,Φ) as a function of B and find the minimum             by inspection.         -   ii. Use a computer program to calculate the loss function             for a range of potential conductor thicknesses, and select             the conductor thickness which yields minimum or desired loss             function value.         -   iii. For the case of R=0, B_(optimum)=π/2.         -   iv. For the case of R=−1 and Φ=0, B_(optimum)=π.         -   v. For the case of R=1 and Φ=0, B_(optimum)=0 (no             incremental macroscopic current in the conductor layer).         -   vi. For the case of R>1 and Φ=0, determine B using the             equation: cos h(B)=R*cos(B).         -   vii. For the case of R>1 and cos(Φ)>0, determine B using the             equation: [R²+1]/R=[cos(B)/cos h(B)+cos h(B)/cos(B)]*cos(Φ).         -   viii. For the general case, determine B using the function:             R cos(Φ)sin(B)[sin h(B)[4 cos²(B)+1]+sin             h(3B)]−[R²+1]sin(2B) sin h(2B)=0. B_(optimum) is given by             the first positive zero of this function. In each case,             physical and/or manufacturing constraints, among other             things, may affect the selection of B, and, therefore,             winding layer thickness.

FIG. 4 is a chart representing an example of the application of this embodiment to determine winding layer thickness in a magnetic component having three windings. More specifically, this example relates to a transformer having one primary (Winding 3) and two secondary (Winding 1 and Winding 2) windings, with an assumed load current in each secondary winding of 1 Amp; an assumed primary excitation current at no load of 1 Amp; a winding length in each winding of 0.1 m (corresponding to element 1 in FIG. 1); six winding layers; and six turns per winding layer. Current phase in Windings 1, 2 and 3 is −1.3, 0.5 and 2.23 radians, respectively, with the phase in Winding 3 derived from assumed values in Windings 1 and 2 in consideration of the no load primary excitation current. From the secondary load currents, the magnitude of the current in the primary winding (Winding 3), 1.88 Amps, is calculated using the condition of load Ampere-Turn equivalence to determine the primary load current component. The total primary current is the vector sum of the primary load current component and the primary no load excitation current component. B, the ratio between winding layer thickness and skin depth, is initially assumed to be 2. The magnitude and phase of peak magnetic surface field intensities below and above each layer (H_(below) and H_(above) in FIG. 4) are determined from the following expressions:

${\overset{\rightarrow}{H}}_{n} = {{\overset{\rightarrow}{H}}_{n - 1} - {\frac{N_{n}}{l}{\overset{\rightarrow}{I}}_{n}}}$ ${\overset{\rightarrow}{H}}_{0} = \frac{\sum\limits_{i = 1}^{n}\;{N_{i}{\overset{\rightarrow}{I}}_{i}}}{2l}$ wherein {right arrow over (H)}_(n) equals H_(n) cos(ωt+φ_(n)) and {right arrow over (I)}_(n) equals I_(n) cos(ωt+Θ_(n)), ω is the radial frequency and φ_(n) and Θ_(n) are phase values. FIG. 2 illustrates {right arrow over (H)}_(n) for an example winding layer 10. Notably, in FIG. 4, the magnitude and phase of the magnetic surface field intensity above Layer Number 1, for example, is equal to the intensity below Layer Number 2, as expected.

The ratio of opposing magnetic surface field intensities (H_(above)/H_(below) in the example) is then computed as well as phase shift. In the example of FIG. 4, normalized dissipation values are then computed for each winding layer. These normalized values are determined from the power loss function, f(H,R,B,Φ), defined in FIG. 4. The dissipation of each element is referenced to a common value. Therefore, in this example, the normalized dissipation has a value of 100%.

To determine the desired winding layer configuration, the expression in step e.vii. above, for the general case, was plotted. The improved (desired or optimum) value for B was then determined as the value for B which caused the first positive zero of the expression. The bottom portion of the chart in FIG. 4 sets forth the winding layer thicknesses (in terms of B) and related boundary conditions for the improved magnetic component. Notably, each winding layer has a different thickness and the power dissipation is reduced by 16.7 percent. In this example, number of layers and turns per layer were held fixed so boundary conditions are unchanged. Further, no specific skin depth was presented since the utilized expression is not dependent upon skin depth, only B. Although true power is indeed a function of skin depth, actual skin depth is irrelevant because the results in this example are displayed in a normalized comparison.

In another preferred embodiment, corresponding to nonsinusoidal excitation having Fourier harmonic components (which may include DC) and fixed winding layers and a fixed number of turns per winding layer, one or more of the following steps are performed to determine a desired winding layer thickness:

-   -   a. Identify the current waveform for each winding in the         magnetic component.     -   b. Using Fourier decomposition, evaluate the harmonic components         of each current waveform, noting magnitude and phase for each         harmonic component.     -   c. Using computer iteration, vary the thickness of each winding         layer to determine the conductor thickness which minimizes total         harmonic dissipation in each layer.         -   i. Using the respective harmonic components of each winding             current, apply Ampere's Law to evaluate the peak magnetic             surface field intensity and phase at each winding layer             boundary for each harmonic frequency.         -   ii. For each layer, and at each harmonic frequency,             determine the ratio of opposing magnetic surface field             magnitudes and the relative phase displacement.         -   iii. For each layer, calculate the respective harmonic             dissipations using the loss function, f(H,R,B,Φ), applied to             each harmonic boundary condition. For a DC component of             current, I_(DC), in a winding layer of N turns and winding             thickness T, the normalized DC power loss function is:

$f_{DC} = {\frac{1}{T}\left( \frac{N + I_{DC}}{l} \right)^{2}}$

-   -   -   iv. Calculate the total harmonic dissipation in each layer             by summing the dissipations for each harmonic frequency.         -   v. For each layer, determine the thickness which minimizes             the total harmonic dissipation.             Each of these steps (as with the steps in the other             disclosed embodiments) may be performed with the aid of a             computer, or through computer software or code, and may be             performed manually.

In another preferred embodiment, corresponding to nonsinusoidal excitation having Fourier harmonic components (including DC) and a variable number of winding layers and a variable number of turns per winding layer, one or more of the following steps are performed to determine desired winding parameters comprising winding layer thickness, the number of winding layers and the number of turns per winding layer:

-   -   a. Identify the current waveform for each winding in the         magnetic component.     -   b. Using Fourier decomposition, evaluate the harmonic components         of each current waveform, noting magnitude and phase for each         harmonic component.     -   c. Using computer iteration, vary the number of layers, the         number of turns per layer, and the conductor thickness of each         winding layer to determine the configuration which minimizes         total winding dissipation.         -   i. For each considered combination of winding layer(s) and             number(s) of turns per layer, apply Ampere's Law to evaluate             the magnetic surface field intensity and phase at each             conductor layer boundary for each harmonic frequency, using             the respective harmonic components of each winding current.         -   ii. For each considered winding layer(s) and number(s) of             turns per layer, and at each harmonic frequency, determine             the ratio of opposing surface field magnitudes and the             relative phase displacement.         -   iii. For each considered combination of winding layer(s) and             number(s) of turns per layer, calculate the respective             harmonic dissipations using the loss function, f(H,R,B,Φ),             applied to each harmonic boundary condition. For a DC             component of current, the normalized power loss function is             expressed as f_(DC) above.         -   iv. For each considered combination of layer(s) and             number(s) of turns per layer, calculate the total harmonic             dissipation by summing the dissipations for each harmonic             frequency.         -   v. For each considered combination of layer(s) and number(s)             of turns per layer, determine the thickness which minimizes             the total harmonic dissipation.         -   vi. Evaluate the minimum dissipation for all other             considered combinations of winding layers and number of             turns per layer using this method.         -   vii. Determine the particular winding configuration (number             of layers, respective turns per layer, and respective             conductor layer thicknesses) which yields minimum total             loss. As in each embodiment, the winding parameters selected             may be limited by cost, manufacturing or physical             constraints, or by a specification such as leakage             inductance or capacitance.

FIG. 5 is a chart representing an example of the application of this embodiment to determine a desired number of winding layers, turns per layer and winding layer thickness for a transformer having one primary (Winding 2) and one secondary (Winding 1) winding. This problem assumed a variable number of winding layers and turns per winding layer. The example also assumed a load current of 1 Amp and a corresponding current in the primary winding of 1.41 Amps considering an assumed primary no load current of 1 Amp. Winding length l is 0.1 meters. The current phase in Windings 1 and 2 is 0.0 and 2.36 radians, respectively. In this example, the method of this embodiment was applied to only one winding (Winding 2). B was initially assumed to be 1.57 for each of two winding layers. The desired configuration was determined as follows: A specific configuration of number of layers and turns per layer was selected. For the resultant corresponding boundary conditions of R and Φ, the desired winding layer thickness was determined by plotting the expression in step e.vii. The improved (desired or optimum) value for B was determined as the value for B which caused the first positive zero of the expression.

The bottom portion of the chart in FIG. 5, labeled “After Improvement,” sets forth the winding layer parameters and related boundary conditions for the improved magnetic component, after applying the method of this embodiment. In this example, power dissipation is improved by 55 percent. One significant aspect of the resulting improved configuration is that the number of turns per layer differs in the three resulting layers, and dissipation is significantly improved. This alone distinguishes from the prior art.

While there has been described what is believed to be the preferred embodiment of the present invention, those skilled in the art will recognize that other and further changes and modifications may be made thereto without departing from the spirit of the invention. For example, the method of the present invention may be applied irrespective of conductor geometry and manufacturing or other physical or cost constraints. Further, all or portions of the inventive method may be applied to all or portions of a magnetic component. Therefore, the invention is not limited to the specific details and representative embodiments shown and described herein. Accordingly, persons skilled in this art will readily appreciate that various additional changes and modifications may be made without departing from the spirit or scope of the invention, as defined and differentiated by the following claims. In addition, the terminology and phraseology used herein is for purposes of description and should not be regarded as limiting. 

1. A method of calculating desired winding parameters for use in a magnetic component having a winding current, comprising: determining the harmonic components of a winding current; iteratively generating a plurality of combinations of one or more of the following winding parameters: number of winding layers, number of turns per winding layer, and winding layer thickness; for one or more of said harmonic components, computing power dissipation for a plurality of said combinations as a function of the relative phase displacement between the magnetic field at the inner boundary and the corresponding outer boundary of one or more of said iteratively generated winding layers; comparing said power dissipations; and selecting a desired combination of winding parameters having desired power dissipation.
 2. The method as set forth in claim 1 wherein said power dissipation is determined as a function of R, Φ and conductor thickness, wherein R is a ratio between the magnetic field intensities between corresponding inner and outer boundaries, and Φ is said relative phase displacement.
 3. The method as set forth in claim 2 wherein said power dissipation is computed in accordance with the power loss function, f(H,R,B,Φ).
 4. The method as set forth in claim 1, further comprising: for one or more of said combinations, plotting the power loss function, f(H,R,B,Φ), for one or more of said harmonic components.
 5. The method as set forth in claim 1 wherein at least one of said iteratively generated winding parameters falls within a range corresponding to at least one manufacturing constraint on at least one physical characteristic of a magnetic component.
 6. The method as set forth in claim 1 wherein said harmonic components comprise a DC component.
 7. The method as set forth in claim 1, further comprising: manufacturing a magnetic component having said desired winding parameters.
 8. A magnetic component made according to the method of claim
 1. 9. The magnetic component as set forth in claim 8 wherein said magnetic component is selected from the group consisting of: coils, inductors, transformers and motors.
 10. A method of calculating desired winding parameters for use in a magnetic component having a winding current, comprising: determining the harmonic components of a winding current; iteratively generating a plurality of combinations of one or more of the following winding parameters: number of winding layers, number of turns per winding layer, and winding layer thickness, for one or more of said harmonic components, computing power dissipation for a plurality of said combinations as a function of the ratio between the magnetic field intensity at the inner boundary and the corresponding outer boundary of one or more of said iteratively generated winding layers; comparing said power dissipations; and selecting a desired combination of winding parameters having desired power dissipation.
 11. The method as set forth in claim 10 wherein said harmonic dissipation is determined as a function of R and Φ, wherein R is said ratio between magnetic field intensities, and Φ is the relative phase displacement between the magnetic field at the inner boundary and the corresponding outer boundary of said iteratively generated winding layers.
 12. The method as set forth in claim 10 wherein said harmonic dissipation is determined as a function of R, Φ and conductor thickness, wherein R is said ratio between magnetic field intensities, and Φ is the relative phase displacement between the magnetic field at the inner boundary and the corresponding outer boundary of said iteratively generated winding layers.
 13. The method as set forth in claim 12 wherein said power dissipation is computed in accordance with the power loss function, f(H,R,B,Φ).
 14. The method as set forth in claim 10, further comprising: for one or more of said combinations, plotting the power loss function, f(H,R,B,Φ), for one or more of said harmonic components.
 15. The method as set forth in claim 10 wherein at least one of said iteratively generated winding parameters falls within a range corresponding to at least one manufacturing constraint on at least one physical characteristic of a magnetic component.
 16. The method as set forth in claim 10 wherein said harmonic components comprise a DC component.
 17. The method as set forth in claim 10, further comprising: manufacturing a magnetic component having said desired thickness.
 18. A magnetic component made according to the method of claim
 10. 19. The magnetic component as set forth in claim 18 wherein said magnetic component is selected from the group consisting of: coils, inductors, transformers and motors.
 20. A method of calculating desired winding parameters for use in a magnetic component having a winding current, comprising: determining the harmonic components of a winding current; iteratively generating a plurality of combinations of one or more of the following winding parameters: number of winding layers, number of turns per winding layer, and winding layer thickness; for one or more of said harmonic components: determining the magnetic field intensity at the inner boundary of one or more of said iteratively generated winding layers in a plurality of said iteratively generated combinations; determining the magnetic field phase at the inner boundary of one or more of said iteratively generated winding layers in a plurality of said iteratively generated combinations; determining the magnetic field intensity at the outer boundary of one or more of said iteratively generated winding layers in a plurality of said iteratively generated combinations; determining the magnetic field phase at the outer boundary of one or more of said iteratively generated winding layers in a plurality of said iteratively generated combinations; calculating ratios between corresponding magnetic field intensities; calculating relative phase displacements between corresponding phases; and determining a desired combination of winding parameters as a function of said calculated ratios and said relative phase displacements.
 21. The method as set forth in claim 20, further comprising: determining the harmonic dissipation in one or more of said iteratively generated winding layers in a plurality of said iteratively generated combinations.
 22. The method as set forth in claim 21 wherein said harmonic dissipation is determined as a function of said ratios between corresponding magnetic field intensities.
 23. The method as set forth in claim 21 wherein said harmonic dissipation is determined as a function of said relative phase displacements between corresponding phases.
 24. The method as set forth in claim 21, further comprising: selecting a desired combination of winding parameters having a desired harmonic dissipation.
 25. The method as set forth in claim 21 wherein said harmonic dissipation is determined as a function of R, Φ and conductor thickness, wherein R is a ratio between the magnetic field intensities between corresponding inner and outer boundaries, and Φ is said relative phase displacement.
 26. The method as set forth in claim 25 wherein said power dissipation is computed in accordance with the power loss function, f(H,R,B,Φ).
 27. The method as set forth in claim 20, further comprising: for one or more of said iteratively generated combinations, plotting the power loss function, f(H,R,B,Φ), for one or more of said harmonic components.
 28. The method as set forth in claim 20 wherein said harmonic components comprise a DC component.
 29. The method as set forth in claim 20, further comprising: manufacturing a magnetic component having said desired winding parameters.
 30. A magnetic component made according to the method of claim
 20. 31. The magnetic component as set forth in claim 30 wherein said magnetic component is selected from the group consisting of: coils, inductors, transformers and motors.
 32. A method of calculating desired winding parameters for use in a magnetic component having one or more winding layers, comprising: varying one or more of the following winding parameters: number of winding layers, number of turns per winding layer, and winding layer thickness; determining boundary conditions at the inner boundary and the outer boundary of one or more winding layers for a plurality of combinations of said varied winding parameters; comparing said boundary conditions at said inner boundaries and said outer boundaries; computing power dissipation as a function of the relationships between said compared boundary conditions.
 33. The method as set forth in claim 32, further comprising: selecting a desired number of winding layers having a desired power dissipation.
 34. The method as set forth in claim 32, further comprising: selecting a desired number of turns per winding layer having a desired power dissipation.
 35. The method as set forth in claim 32, further comprising: selecting a desired winding layer thickness having a desired power dissipation.
 36. The method as set forth in claim 32 wherein said power dissipation is computed in accordance with the power loss function, f(H,R,B,Φ).
 37. The method as set forth in claim 32 wherein said comparing comprises plotting the power loss function, f(H,R,B,Φ), for a plurality of said combinations.
 38. The method as set forth in claim 32, further comprising: manufacturing a magnetic component having said desired winding parameters.
 39. A magnetic component made according to the method of claim
 32. 40. The magnetic component as set forth in claim 39 wherein said magnetic component is selected from the group consisting of: coils, inductors, transformers and motors.
 41. A system for calculating desired winding parameters for use in a magnetic component, said system determining power dissipation as a function of one or more boundary conditions at the inner boundaries and the outer boundaries of the winding layers of said magnetic component, and employing the steps of determining desired winding parameters comprising at least one selected from the group consisting of: number of winding layers, number of turns per winding layer, and winding layer thickness.
 42. The system as set forth in claim 41 wherein said boundary conditions comprise at least one of the following: magnetic field intensity and magnetic field phase.
 43. The system as set forth in claim 41, further comprising: selecting a desired number of winding layers having a desired power dissipation.
 44. The system as set forth in claim 41, further comprising: selecting a desired number of turns per winding layer having a desired power dissipation.
 45. The system as set forth in claim 41, further comprising: selecting a desired winding layer thickness having a desired power dissipation.
 46. The system as set forth in claim 41 wherein said power dissipation is determined in accordance with the power loss function, f(H,R,B,Φ).
 47. The system as set forth in claim 41 wherein the power loss function, f(H,R,B,Φ), is plotted for a plurality of combinations of said winding parameters.
 48. A magnetic component made according to the method of claim
 41. 